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SECTION  1 


INTRODUCTIOJ 


This  Study  establishes  guidelines  for  determining  how  accurately  aircraft 
position  and  velocity  should  be  measured  on  a  Tactical  Aircrew  Combat  Training 
System  (TACTS)  range.  The  method  of  analysis  is  based  on  the  sensitivities  of 
Measures  of  Effectiveness  (MOE)  computed  from  TACTS  weapon  simulations  to  input 
aircraft  position  and  velocity  at  the  time  of  weapon  launch.  On  the  basis  of 
desired  MOE  accuracies,  maximum  allowable  errors  (standard  deviations)  are 
determined  for  launch  aircraft  position  and  velocity  estimates  for  a  variety  of 
weapoiVtarget  scenarios. 
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SECnOJ  2 


APPROACH 


TACTS  ranges  support  the  training  of  Navy  pilots.  One  range  function  is  to 
track  the  status  of  each  aircraft,  including  position,  speed,  and  attitude. 
The  accuracy  in  estimated  values  for  these  parameters  depends  upon  many 
factors,  including  geography,  hardware,  telemetry,  and  software  for  data 
processing. 

What  accuracies  are  required?  In  part,  the  answer  lies  in  subsequent  uses 
of  parameter  estimates.  One  such  use  is  an  input  to  weapon  simulation  programs 
;diich  simulate  weapon  performance  and  calculate  appropriate  MOEs. 

Aircraft  position  is  typically  specified  by  values  for  three  coordinates, 
X,  y,  Z,  and  speed  is  specified  three  velocity  components,  V  ,  V  ,  V  .  The 
following  analysis  focuses  on  the  required  accuracy  of  thes^  petition  and 
velocity  conponent  values. 

Suppose  there  were  no  errors  in  the  determination  of  position  and  velocity. 
Then,  assuming  that  the  weapon  simulation  programs  are  completely  accurate,  a 
correct  NOE  value  would  be  computed  in  every  instance.  In  actuality,  there  are 
errors  in  position  and  velocity  estimates,  and  consequently  weapon  simulations 
do  not  always  conpute  an  accurate  MOE  value.  In  general,  the  larger  the  errors 
in  input  position  and  velocity  values,  the  larger  will  be  the  error  in  computed 
MOE.  The  actual  error  sizes  are  not  known  in  advance;  if  so,  they  could  be 
eliminated  in  the  processing  of  data.  At  best,  error  sizes  can  be  character¬ 
ized  by  using  the  language  of  probability  (i.e.,  small  errors  inply  highly 
accurate  MOEs). 

To  formulate  the  analytical  method  for  the  determination  of  TACTS  range 
accuracy  requirements  the  following  guidelines  were  established; 

o  The  method  should  be  general  enough  to  apply  to  various  aircraft/ 
weapon/target  combinations. 

o  The  method  should  require  only  readily  available  input. 

o  The  method  should  produce  quantitative  output  which  would  support  trade¬ 
off  eunalysis  and  provide  flexibility  to  the  decision-^king  process. 

The  goal  of  this  method  was  to  estimate  parametrically  the  amounts  by  which 
aircraft  position  and  velocity  component  values,  as  determined  by  TACTS  range 
telemet^  and  processing,  could  be  in  error  \diile  requiring  that  TACTS  weapon 
simulations  conpute  a  reasoneibly  accurate  MOE  value  reasonably  often. 

As  a  first  step  in  developing  the  analytical  method,  a  Performance  Goal  was 
defined  as  follows: 

Performance  Goal  -  The  accuracy  with  which  the  MOE  is  to  be  determined. 
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Maxinrum  allowable  uncertainties  in  position  and  velocity  component  values 
consistent  with  achievement  of  this  Performance  Goal  were  determined.  For 
example,  for  target-point  target  and  weapon-bomb,  probability  of  kill  (P.  )  is 
an  MOE.  A  Performance  Goal  might  be  a  relative  error  of  no  more  than  ten 
percent  (i.e.,  the  difference  between  confuted  and  correct  values  for  P.  should 
be  no  larger  than  ten  percent  of  the  correct  value  for  P.  ).  Achievement  of 
this  goal  places  restrictions  on  the  magnitudes  of  errors  in  the  values  of 
position  and  velocity  components. 

The  errors  in  position  and  velocity  components  are  treated  as  unbiased, 
independent,  normal  random  variables.  Qualitatively  this  means  that  errors  are 
just  as  likely  to  be  positive  as  negative  with  the  average  error  being  zero. 
Further,  the  error  in  determining  one  position  or  velocity  component  value  is 
not  related  to  the  error  in  determining  a  different  conponent. 

For  errors  v^ich  are  normal,  the  likelihood  of  errors  is  determined  by  the 
distribution  parameter  "standard  deviation"  or  sigma  ( a) .  Hie  larger  the  value 
of  a,  the  more  likely  will  be  the  occurrence  of  large  errors.  The  magnitudes 
of  position  and  velocity  component  errors  for  achievement  of  a  Performance  Goal 
can  be  stated  in  terms  of  error  distribution  standard  deviations. 

A  requirement  that  the  Performance  Goal  must  be  achieved  almost  always  will 
require  that  TACTS  range  position  and  velocity  estimates  be  very  accurate 
al^st  always.  On  the  other  hand,  if  achieving  the  Performance  Goal  less 
frequently  is  acceptable,  then  TACTS  range  estimates  must  be  accurate  less 
frequently  (and  consequently  accuracy  requirements  would  be  relaxed).  Thus,  we 
define  a  Performance  Level  as  follows: 

Performance  Level  -  The  frequency  (fraction  of  launches)  that  the  Perform- 

ance  Goal  is  achieved. 

Performance  Levels  are  varied  from  .75  to  .95  in  this  analysis.  Both  GDOP 
<5  and  5<  ODOP  <12  are  considered  (QDOP  stands  for  Geometric  Dilution  of 
Precision,  a  function  of  the  TACTS  range  physical  layout  of  position-locating 
equipment,  reference  a).  Under  favorable  conditions  (QX)P  <5),  vertical 
position  U)  and  velocity  (V)  components  are  assumed  to  be  determined  with  the 
same  accuracy  as  horizontaf  (V'^  and  Y/V  )  values.  under  less  favorable 
conditions  (5<  (HXDP  <12 ),  standard* deviations  for  errors  in  vertical  component 
determinations  are  assumed  to  be  four  times  as  large  as  the  standard  deviations 
in  horizontal  conponent  determinations. 

For  each  combination  of  target,  weapon,  and  NOE,  a  parametric  set  of  TACTS 
range  accuracy  requirements  is  generated  by  the  following  general  procedure*: 

1.  Specify  the  Performance  Goal:  the  accuracy  with  which  the  MOE  is  to  be 
determined  (e.g.,  computed  Pj^  must  be  within  ten  percent  of  act\ial  Pj^). 

2.  Identify  values  for  position  and  velocity  components  reconmended  by 
tactical  considerations  (these  are  "intended"  values  at  weapon 
release ) . 


*  E)ocinnented  more  fully  in  .^pendix  A. 
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3. 


Using  the  weapon  simulation: 

a.  Simulate  weapon  lamch  with  "intended"  values  and  move  target  to 
weapon  inqpact  position. 

b.  For  each  of  several  variations**  (representing  "pilot  error")  in 
position  and  velocity  conqponent  values  about  their  "intended" 
values : 

(1)  Simulate  weapon  launch,  determining  conpjted  MOE  (the  "cor¬ 
rect"  value). 

(2)  In  each  of  the  six  components  (three  position  and  three 
velocity),  one  at  a  time,  determine  the  extent  to  vAiich  the 
value  can  be  changed  while  still  producing  a  computed  MOE 
value  vAiich  is  within  the  accuracy  of  the  "correct"  value, 
specified  in  the  Performance  Goal  (Step  1).  This  determina¬ 
tion  is  done  by  running  the  weapon  simulation  with  altered 
input  values. 

c.  At  this  stage,  several  conponent  value  intervals  have  been  deter¬ 
mined  vrfiich  represent  allowable  variability  xander  the  accuracy 
requirement  specified  in  Step  1. 

d.  Repeat  from  Step  2  if  there  are  additional  tactical  cases  of 
importance. 

4.  Specify  the  Performance  Level  (i.e.,  fraction  of  time  that  the  computed 
MOE  value  is  to  be  within  the  accuracy  specified  for  the  Performance 
Goal  in  Step  1). 

5.  Using  the  BASIC  program  DISCRETE  FIT***,  compute  combinations  of 
standard  deviations  (<j's)  for  position  and  velocity  component  values 
for  which  errors  within  the  bounds  identified  in  Step  3.c  occur  with 
the  frequency  (fraction  of  time)  specified  in  Step  4.  The  resulting 
<y's  constitute  TACTS  range  parametric  accura^  requirements  produced  by 
this  methodology  for  the  selected  combination  of  target,  weapon,  and 
MOE  at  the  specified  Performance  Level. 

6.  Repeat  from  Step  4  for  additional  Performance  Levels. 

To  summarize,  the  TACTS  weapon  simulations  were  used  to  determine  in 
absolute  terms  how  much  variability  can  be  tolerated  in  position  and  velocity 
conponent  determination  v4iile  still  maintaining  acceptable  MOE  ac^racy. 
Treating  these  absolute  quantities  as  random  variables,  the  statistical 
approach  was  to  conpate  corresponding  standard  deviations  v4iich  assure  that 
tolerable  errors  in  position  and  velocity  conqponent  determination  are  achieved 
with  specified  frequency. 

Presented  in  the  following  sections  is  an  extensive  collection  of  computed 
results  pertaining  to  the  following  weapon/target  combinations: 


**  See  Appendix  B. 

***  See  Appendix  C. 
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1.  Bombing  a  point  target  (three  cases) 

2.  Bombing  a  runway  (one  case) 

3.  Bombing  a  point  target  -  retarded  bomb  (one  case) 

4.  Firing  a  rocket  at  a  point  target  (three  cases) 

5.  Firing  a  gun  in  air-to-air  combat  (seven  cases) 

The  different  cases  considered  in  each  scenario  represent  different  intended 
weapon  launch  conditions  reflecting  tactical  options.  Cases  are  defined,  case 
data  is  reported,  scenario  MOEs  and  Performance  Goals  are  given,  and  confuted 
results  are  presented  graphically,  parametrized  by  Performance  Level.  Illus¬ 
trations  of  the  interpretation  and  use  of  these  parametrically  presented 
resxilts  are  provided. 
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SECTION  3 


TACTS  RANGE  ACCURACY  GUIDELINES 
Bombing  a  Point  Target  Scenario 


The  following  conditions  are  assumed  for  the  point  target  bombing  scenario: 

Target  -  Small  fighter  aircraft  on  the  ground 

Weapon  -  MK-82  Mod  1  (Std)  electric  fuse  bomb 

Aircraft  -  F/A-18 

MOE  -  Pj^ 

Performance  Goal  -  No  more  than  ten  percent  relative  error  (alternately,  no 

more  than  five  percent  relative  error). 

Performance  Levels  -  .75,  .80,  .85,  .90,  .95 

Presented  in  Table  1  are  input  values  for  position,  speed,  pitch,  heading, 
and  the  velocity  components  conpjted  from  the  latter  dynamic  parameters  for  the 
three  cases  investigated  in  t^is  scenario.  The  bomb  drop  parameters  were 
chosen  from  a  Tactical  Manual  (reference  b). 

In  order  to  simulate  pilot  error  in  the  attainment  of  "intended"  launch 
parameter  values,  the  method  documented  in  Appendix  B  was  employed  which 
yielded  the  values  presented  in  Table  1  for  each  case  investigated. 

Figures  1  and  2  present  in  graphical  and  parametric  form  scenario-specific 
guidance  for  assessing  current  TACTS  range  accuracies  or  estedslishing  future 
accuracy  requirements.  Figures  1  and  2  pertain  to  (aX)P  <5  and  to  5<  GDOP  <12, 
res^ctively.  Figures  3  and  4  present  guidance  derived  under  the  more 
stringent  Performance  Goal  that  there  be  no  more  than  five  percent  relative 
error  in  the  determination  of  Pj^. 

To  use  the  figures,  first  select  a  value  for  V,  the  standard  deviation  for 
velocity  component  error,  on  the  horizontal  axis,  and  a  value  for  L,  the 
standard  deviation  for  position  component  error,  on  the  vertical  axis.  The 
corresponding  point  (V,L)  defines  the  Performance  Level;  for  instance,  if  (V,L) 
falls  within  (below)  the  curve  ledselled  .85,  then  for  TACTS  reuige  accuracies  L 
and  V,  the  Performance  Level  will  be  at  least  85  percent.  If  (V,L)  falls 
outside  (above)  the  curve  labelled  .85,  then  for  TACTS  range  accuracies  L  and 
V,  the  Performance  Level  will  be  less  than  85  percent.  Further,  all  (V,L) 
pairs  within  the  .85  contour  will  yield  Performance  Levels  of  at  least  85 
percent. 

Figures  1  through  4  provide  a  ready  vehicle  for  establishing  required  TACTS 
reinge  accuracies  for  a  single  target  type  under  <3X)P  <5  or  5<  <3X)P  <12  condi¬ 
tions.  For  instance,  consider  the  following  typicar  conclusions  obtainable 
from  Figure  1: 
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1.  If  telemetry  leads  to  a  standard  deviation  in  aircraft  velocity  conpo- 
nents  of  V  -  1  ft/sec  and  a  standard  deviation  in  aircraft  position 
components  of  L  «  5  ft  then  the  bomb  simulation  would  produce  a  value 
for  P.  v^ose  relative  error  is  no  more  than  10  percent  between  85 
percem:  and  90  percent  of  the  time. 

2.  If  velocity  conponent  error  standard  deviations  are  V  >  1  ft/sec,  then 
the  Performance  Level  can  never  be  as  large  as  87.5  percent. 

3.  If  position  component  error  standard  deviations  are  L  ■  1  ft,  then 
velocity  conponent  error  standard  deviations  (V)  must  not  exceed 
approximately  6  ft/sec  if  a  Performance  Level  of  85  percent  is 
required. 
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TABLE  1.  CASE  DATA  -  SIMULATED  BOMB  DROPS  ON  A  POINT  TARGET 


Case  12  3 

(1)  Intended  aircraft  position  relative  to  aim  point  on  target  at  time  of  bomb 
drop 

Xpos  4129  ft  8346  ft  4790  ft 

Ypos  0  ft  0  ft  0  ft 

Zpos  5000  ft  2000  ft  2000  ft 

(2)  Intended  dynamic  parameter  values  at  time  of  bomb  drop 

Speed  802_ft/sec  760^ft/sec  760„ft/sec 

Pitch  -45°  0°  -15° 

Heading  0  0  0° 

(3)  Velocity  con?)onents  computed  from  (fynamic  data 


567.1  ft/sec 

760 

ft/sec 

724.44  ft/sec 

vj 

0  ft/sec 

0 

ft/sec 

0  ft/sec 

-567.1  ft/sec 

0 

ft/sec 

-194.12  ft/sec 

Sigmas  conpated  for  pilot  error 

ax 

17.8  ft 

41.85 

ft 

22.75  ft 

oy 

19.8  ft 

25.9 

ft 

15.81  ft 

oz 

24.3  ft 

20.3 

ft 

12.83  ft 

2.35  ft/sec 

3.87 

ft/sec 

3.48  ft/sec 

2 . 7  ft/sec 

2.35  ft/sec 

2 . 4  ft/sec 

< 

3 . 3  ft/sec 

1.85 

ft/sec 

1.95  ft/sec 
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Notes:  V  -  Velocity  Component  Error  Standard  Deviation  (ft/sec) 

(1)  ~  ^  T 

Fiqure  1.  TACTS  Range  Accuracy  -  Bomb  vs.  Point  Target 

(2)  V-oVj^=ciV  -aVjj,  (GDOP  <  5,  10%  Relative  Error  in  Pj^) 


SECTIOJ  4 


TACTS  RANGE  ACCURACY  GUIDELINES 
Bombing  a  Runway  Scenario 


The  following  conditions  are  assvuned  for  the  ninway  bombing  scenario: 

Target  -  Long  runway 

Weapon  -  MK-82  Mod  1  (Std)  electric  fuse  bomb 

Aircraft  -  P/A-18 

MOE  -  Crater  Diameter 

Performance  Goal  -  Compute  correct  crater  diameter  exactly 

Performance  Levels  -  .75,  .80,  .85,  .90,  .95 

Presented  in  Table  2  are  input  values  for  position,  s^ed,  pitch,  heading, 
and  the  velocity  conponents  conputed  from  the  latter  dynamic  parameters  for  the 
case  investigated  in  this  scenario.  The  parameters  were  chosen  from  a  Tactical 
Manual  (reference  b).  In  addition,  sigmas  computed  for  pilot  error  are  listed 
for  the  position  and  velocity  conponents. 

Figures  5  and  6  present  in  graphical  and  parametric  form  scenario-specific 
guidance  regarding  TACTS  range  accuracy.  Figures  5  and  6  pertain  to  GDOP  <5 
and  to  5<  GDOP  <12,  respectively. 
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TABLE  2.  CASE  DATA  -  SIMULATED  BOMB  DROPS  (X(  A  RUNWAY 


Case  1 

(1)  Intended  aircraft  position  relative  to  aim  point  on  target  at  time  of  bomb 
drop 

Xpos  6206  ft 

Ypos  0  ft 

Zpos  8000  ft 

(2)  Intended  dynamic  parameter  values  at  time  of  bomb  drop 

Speed  844  ft/sec 

Pitch  -45° 

Heading  0° 

(3)  Velocity  con?x)nents  con^juted  from  dynamic  data 

596.8  ft/sec 
0  ft/sec 
596.8  ft/sec 

(4)  Sigmas  computed  for  pilot  error 


vY 


OX 

25.8  ft 

oy 

30,9  ft 

az 

41.1  ft 

2.4  ft/sec 

dv 

V 

2.9  ft/sec 

z 

3.9  ft/sec 
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dV  /4  Figure  6.  TACTS  Range  Accuracy  -  Bomb  vs.  Runway  Target 

^  (5  <  GDOP  <  12,  Correct  Crater  Diameter) 


SECTIC»I  5 


TACTS  RANGE  ACCURACY  GUIDELINES 
Bombing  a  Point  Target  With  Retarded  Bomb  Scenario 


The  following  conditions  are  assumed  for  the  bombing  of  a  point  target  with 
a  retarded  bomb  scenario: 

Target  -  Small  fighter  aircraft  on  the  ground 

Weapon  -  MK~82  Snakeye  retarded  bomb 

Aircraft  -  F/A-18 

MOE  - 

Performance  Goal  -  No  more  than  ten  percent  relative  error 
Performance  Levels  -  .75,  .80,  .85,  .90,  .95 

Preseflted  in  Table  3  are  input  values  for  position,  speed,  pitch,  heading, 
and  the  velocity  conponents  confuted  from  the  latter  dynamic  parameters  for  the 
one  case  investigated  in  this  scenario.  Pilot  error  sigmas  for  position  and 
velocity  conponents  are  also  shown.  Uie  initial  input  parameters  were  chosen 
from  a  Tactical  Manual  (reference  b).  A  retarded  bomb  was  used  in  this 
scenario  for  con5»tibility  with  the  low^altitude,  higlv-speed  "pop-up"  delivery 
technique . 

Figures  7  and  8  present  in  graphical  and  parametric  form  scenario-specific 
guidemce  regarding  TACTS  reunge  accuracy  guidelines.  Figures  7  and  8  pertain  to 
GDOP  <5  and  to  5<  (3XDP  <12,  respectively. 
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TABLE  3.  CASE  DATA  -  SIMULATED  BOMB  DROPS  OJ  A  POINT 
TARGET  -  RETARDED  BOMB 


Case 


(1)  Intended  aircraft  position  relative  to  aim  point  on  target  at  time  of  bomb 
drop 


Xpos 

Ypos 

Zpos 


1468  ft 
0  ft 
500  ft 


(2)  Intended  dynamic  parameter  values  at  time  of  bomb  drop 


Speed 

Pitch 

Heading 


760.ft/sec 

-15? 


(3)  Velocity  conponents  conputed  from  dynamic  data 


V 


724.44  ft/sec 
0  ft/sec 
-194.2  ft/sec 


(4)  Sigmas  computed  for  pilot  error 


OX 

12.9  ft 

ay 

5.14  ft 

az 

3.78  ft 

aV 

< 

5.16  ft/sec 
2.04  ft/sec 

oV* 

z 

1 . 5  ft/sec 

18 


19 


Ul 


20 


SECTICXJ  6 


TACTS  RANGE  ACCURACY  GUIDELINES 
Rocket  Firing  at  a  Point  Target  Scenario 


The  following  conditions  are  assmed  for  the  rocket  at  a  point  target 
scenario : 

Target  -  Small  fighter  aircraft  on  the  ground 

Weapon  -  MK-4/MK-67  WP/M  427  2.75"  Folding  Fin  Aircraft  Rocket  (FFAR) 

Aircraft  -  F/A-18 

MOE  -  Miss  Distance  (In  current  rocket  simulations,  the  post  impact  scoring 
routine  is  not  invoked;  therefore,  Pj^  is  not  calculated.) 

Performance  Goal  -  No  more  than  ten  percent  relative  error 

Performance  Levels  -  .75,  .80,  .85,  .90,  .95 

Presented  in  Table  4  are  input  values  for  position,  speed,  pitch,  heading, 
and  the  velocity  components  computed  from  the  latter  d^amic  parameters  for  the 
three  cases  investigated  in  this  scenario.  Pilot  error  sigmas  for  position  and 
velocity  components  are  also  shown.  The  initial  parameters  were  obtained  from 
Tactical  Manuals  and  represent  a  varied  set  of  realistic  conditions  (reference 
b). 


Figures  9  and  10  present  in  graphical  and  parametric  form  scenario-specific 
guidance  regarding  TACTS  range  accuracy  guidelines.  Figures  9  and  10  pertain 
to  aX)P  <5  and  to  5<  (3XIP  <12,  respectively. 
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TABLE  4.  CASE  DATA  -  SIMULATED  ROCKET  FIRING  AT  A  POINT  TARGET 


Case  123 

(1)  Intended  aircraft  position  relative  to  aim  point  on  target  at  time  of 
rocket  firing 

Xpos  0  ft  0  ft  0  ft 

Ypos  0  ft  0  ft  0  ft 

Zpos  5000  ft  3000  ft  1000  ft 

(2)  Intended  dynamic  parameter  values  at  time  of  rocket  firing 

S^ed  759.9  ft/sec  675.5  ft/sec  665.2^ft/sec 

Pitch  0  0°  0° 

Heading  -30°  -45°  -10° 

(3)  Velocity  coitponents  conpited  from  dynamic  data 


658.1  ft/sec  477.65  ft/sec 
V"  -379.9  ft/sec  -477.65  ft/sec 
Vg  *  0  ft/sec  0  ft/sec 

(4)  Sigmas  conqputed  for  pilot  error 


665.2  ft/sec 
-117.2  ft/sec 
0  ft/sec 


ox 

69.2 

ft 

76.1 

ft 

120.8 

ft 

oy 

68.3 

ft 

56.2 

ft 

38.6 

ft 

oz 

57.8 

ft 

54.1 

ft 

23.9 

ft 

oV 

7.3 

ft/sec 

5.2 

ft/sec 

33.4 

ft 

3.17 

ft/sec 

3.8 

ft/sec 

2.3 

ft/sec 

oV^ 

z 

2.2 

ft/sec 

2.17 

ft/sec 

1.4 

ft/sec 
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SECTiai  7 


TACTS  RANGE  ACCURACY  GUIDELINES 
Gunfire  In  Air-to-Air  Combat  Scenario 


The  following  conditions  are  assumed  for  the  gunfire  in  air-to-air  combat 
scenario: 

Target  -  Small  fighter  aircraft 

Weapon  -  M61A1  20  mm  gun 

Aircraft  -  F-14 

MOE  -  Pj^ 

Performance  Goal  -  No  more  than  ten  percent  relative  error 

Performance  Levels  -  .75,  .80,  .85,  .90,  .95 

Presented  in  Table  5  are  input  values  for  position,  speed,  pitch,  heading, 
and  the  velocity  components  computed  from  the  latter  dynamic  parameters  for  the 
seven  cases  investigated  in  this  scenario.  Pilot  error  sigmas  for  position  and 
velocity  conponents  are  also  shown.  The  parameters  were  chosen  from  a  Tactical 
Manual  (reference  b).  Sigmas  for  x  position  and  x  velocity  parameters  are  not 
presented  since  the  analysis  indicated  total  NOE  insensitivity  to  euiy  change  in 
x-direction  conponents. 

Figures  11  and  12  present  in  graphical  and  parametric  form  scenario- 
specific  guidance  regarding  TACTS  reinge  accuracy  guidelines.  Figures  11  and  12 
pertain  to  C3X)P  <5  and  to  5<  GDOP  <12,  respectively.  The  standard  deviations 
presented  pertain  to  the  differences  between  fighter  and  target  position  compo¬ 
nents,  and  between  fighter  and  target  velocity  conponents. 

For  scenarios  involving  an  airborne  target,  both  launcher  and  target 
aircraft  position  and  velocity  are  estimated.  Under  the  assunption  that 

estimation  errors  are  independent  and  statistically  the  same  for  both  aircraft, 
the  error  standard  deviations  for  both  the  lawcher  and  target  are  each  1//2 
times  the  error  standard  deviation  in  the  estimate  of  the  difference  of  like 
parameter  values.  For  instance,  if  and  X-^-  are  the  errors  in  estimates  of 

launcher  and  target  X-conponent  values  at  time  of  weapon  firing,  then  x  » 

*^Xl^  +  If  ®Xl^  -  ®x  ^  »  *^x^,  then  ^xi^-x  -  2  ®x^  and  "^x  »  1//2  ‘^Xi-^"''. 

T  T  T  T 


Guidance  for  accuracy  requirements  based  vpon  airborne  target  scenarios  is 
presented  in  Figures  11  and  12  in  terms  of  paired  values  (^v^— v  ,  ^x^— x  ).  If 
a  partioilar  ninnerical  pair  (V,L)  assures  a  Performance  Level  of  .90,  then  the 
required  single  aircraft  error  standard  deviation  pair  is  (V//2,  L//2). 
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TABLE  5.  OVSE  DATA  -  SIMULATED  AIR-TO-AIR  GUNFIRE 


Case 

1 

2 

3 

(1)  Intended  aircraft  position  at 

time  of  firing 

Latincher: 

Xpos 

Ypos 

Zpos 

0  ft 

0  ft 

15000  ft 

0  ft 

0  ft 

15000  ft 

0  ft 

0  ft 
15000  ft 

Target: 

J^s 

Ypos 

Zpos 

1500  ft 

0  ft 

15000  ft 

500  ft 

0  ft 

15000  ft 

1500  ft 

0  ft 
15000  ft 

( 2 )  Intended 

dynamic  parameter  values  at  time  of  firing 

Launcher : 

Speed 

Pitch 

Heading 

500  ft/sec 
■2° 

0° 

500  ft/sec 
"2° 

0° 

500  ft/sec 

6° 

Target: 

Speed 

Pitch 

Heading 

500_ft/sec 

-2° 

0° 

500_ft/sec 

-2° 

500  ft/sec 
0° 

45° 

(3)  Velocity 

components  conpjted  from  dynamic  data 

Launcher: 

"Sc 

499.69  ft/sec 

0  ft/sec 
-17.44  ft/sec 

499.69  ft/sec 

0  ft/sec 
.  -17.44  ft/sec 

497.26  ft/sec 
52.27  ft/sec 

0  ft/sec 

Target : 

V 

z 

500  ft/sec 

0  ft/sec 
-17.44  ft/sec 

500  ft/sec 

0  ft/sec 
-17.44  ft/sec 

353.33  ft/sec 
353.55  ft/sec 

0  ft/sec 

(4)  Sigmas  confuted  for  pilot  error 

oy 

CTZ 

oV 

z 

5.5  ft 

5.6  ft 

12.08  ft/sec 

12.07  ft/sec 

1.84  ft 

1.84  ft 

12.08  ft/sec 

12.07  ft/sec 

5.35  ft 

5.36  ft 

12.09  ft/sec 
12.07  ft/sec 
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TABLE  5.  CASE  DATA  -  SIMULATED  AIR-TO-AIR  GUNFIRE  (Cont.) 


Case  i  5  6  7 

(1)  Intended  aircraft  position  at  time  of  firing 
Launcher: 


}^s 

0  ft 

0  ft 

0  ft 

0 

ft 

Ypos 

0  ft 

0  ft 

0  ft 

0 

ft 

Zpos 

15000  ft 

15000  ft 

15000  ft 

15000 

ft 

Target: 

Xpos 

500  ft 

1500  ft 

500  ft 

1500 

ft 

Ypos 

0  ft 

0  ft 

0  ft 

0 

ft 

Zpos 

15000  ft 

15000  ft 

15000  ft 

15000 

ft 

(2)  Intended  dynamic  parameter  values  at  time  of  firing 
Launcher: 

Speed  500  ft/sec  500„ft/sec  500^ft/sec  500  ft/sec 

Pitch  0°  0°  0°  0° 

Heading  6®  9®  9®  0® 

Target: 

Speed  500„ft/sec  500„ft/sec  500^ft/sec  500  ft/sec 

Pitch  0®  0®  0®  0° 

Heading  45®  90®  90®  180® 

(3)  Velocity  components  confuted  from  dynamic  data 
Launcher: 


V  497.26  ft/sec  493.84  ft/sec 

Vj  52.26  ft/sec  78.21  ft/sec 

V'  0  ft/sec  0  ft/sec 

Target: 

V  353.55  ft/sec  0  ft/sec 

V"  353.55  ft/sec  500  ft/sec 

0  ft/sec  0  ft/sec 

(4)  Sigmas  computed  for  pilot  error 


oy 

1.76 

ft 

4.82 

ft 

az 

1.77 

ft 

4.83 

ft 

12.08 

ft/sec 

12.09 

ft/sec 

(jV^ 

Z 

12.07 

ft/sec 

12.07 

ft/sec 

493.84  ft/sec 
78.21  ft/sec 
0  ft/sec 


0  ft/sec 
500  ft/sec 
0  ft/sec 


1.6  ft 
1.6  ft 
12.08  ft/sec 
12.07  ft/sec 


500  ft/sec 
0  ft/sec 
0  ft/sec 


-500  ft/sec 
0  ft/sec 
0  ft/sec 


4.17  ft 
4.17  ft 
12.08  ft/sec 
12.07  ft/sec 
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SECTIMJ  8 


COMPARISOJS 


TACTS  range  accuracy  guidelines  for  five  different  weapon/target  scenarios 
were  presented  in  the  previous  five  sections.  For  coitparison,  the  90  percent 
Performance  Level  ^rves  under  GDOP  <5  conditions  for  each  of  these  scenarios 
are  presented  in  Figure  13.  Note  that  the  V  and  L  scales  are  NOT  the  same  for 
all  five  curves.  - 

The  most  demanding  scenarios,  in  terms  of  TACTS  range  accuracy  guidelines, 
are  firing  a  gun  in  air— to— air  combat  and  bombing  a  point  target  with  a  re— 
tarded  bomb.  The  least  dem^ding  is  bombing  a  runway.  in  the  former 
scenarios,  small  errors  in  position  and  velocity  conf)onent  values  cause  large 
variations  in  R  .  In  the  latter  scenario,  the  simulation  is  more  tolerant  of 
positron  ^d  velocity  conponent  values.  In  order  of  simulated  weapon 
effectiveness  sensitivity  to  TACTS  range  accuracy,  the  five  scenarios  are 
ranked  from  most  demanding  to  least  demanding  as  follows; 

Firing  a  gun  in  air-to-air  combat  or 
Bombing  a  point  target  with  a  retarded  bomb 
Bombing  a  point  target 
Firing  a  rocket  at  a  point  target 
Branbing  a  runway 

Generally,  the  more  sensitive  weapon  effectiveness  is  in  position  and 
velocity  component  values  at  weapon  launch,  the  more  stringent  are  the  TACTS 
range  accuracy  requirements.  Thus,  the  results  of  this  analysis  do  conform 
qualitatively  with  operational  realities. 

The  TACTS  range  accuracy  requirements  should  NOT  be  combined  by  simple 
averaging.  Any  requirement  consolidation  should  “Include,  at  a  minimum,  a 
detennination  of  the  relative  frequency  with  which  each  scenario  is  to  be 
exercised.  If  the  TACTS  range  were  instrumented  to  conform  with  a  "consoli¬ 
dated  accuracy  requirement,"  then  high  Performance  Levels  would  be  realized  for 
less  demanding  scenarios  and  low  Performance  Levels  for  the  more  demanding 
scenarios.  If  different  scenarios  are  to  be  exercised  on  the  same  portion  of 
the  TACTS  range,  then  requirement  consolidation  is  unavoidable. 

The  results  of  this  stu(fy  also  ceui  be  used  to  assess  how  well  a  particular 
range  instrumentation  would  perform.  For  instance,  suppose  that  TACTS  range 
instrumentation  allows  for  determination  of  position  components  with  error 
st^dard  deviation  of  1  ft  and  velocity  conponents  with  error  standard  devi¬ 
ation  of  2.5  ft/sec.  Then  confuted  values  of  P.  (with  ODOP  <  5)  will  have  a 
relative  error  of  no  more  than  10  percent  for:  (a)  between  75  percent  and  80 
percent  of  all  retarded  bomb  drops  on  point  targets  (Figure  7);  (b)  between  85 
percent  and  90  percent  of  all  standard  bomb  drops  on  a  point  target  (Figure  1); 
and  (c)  between  75  percent  and  80  percent  of  all  gun  firings  against  an 
airborne  target  (Figure  11).  ^  ^ 
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Gunfire/Aircraft  2.6  +  Retarded  Bomb/Point  Target  I  Bomb/Point  Target 


jure  13.  TACTS  R2«ige  Accuracy  Requirements  Determined  frcan  Five  Weapon/Target  Combinations 

(90%  Performance  Level/5  <  GEXDP) 


Some  options  for  determining  specific  TACTS  range  accuracy  requirements  are 
as  follows: 

a.  Accommodate  the  most  stringent  accuracy  requirements  to  satisfy  a 
desired  Performance  Level  for  the  most  demanding  scenario.  Then 
desired  Performance  Levels  for  other  scenarios  will  be  exceeded. 

b.  Irtplement  the  best  accuracy  requirements  consistent  with  technological 
and  financial  resources.  Then  there  could  be  a  disparity  in  Perfor¬ 
mance  Levels  for  certain  scenarios. 

c.  Implement  different  accuracy  requirements  for  different  multi-  and 
single-purpose  sites  on  the  TACTS  range.  Then  range-wide  complicince 
with  Performance  Level  requirements  might  be  achievable  within  the 
constraints  of  technological  and  financial  resources. 
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SECTIOJ  9 


COICLUSIWS 


Specification  of  TACTS  range  accuracy  will  depend  in  part  upon  available 
range  instrumentation  technology.  The  figures  presented  provide  a  vehicle  for 
trading  off  position  and  velocity  accuracy  requirements  for  areas  of  the  range 
used  for  different  weapon/target  pairs. 

Large  values  for  Performance  Level  (.90  to  .95)  lead  to  stringent  TACTS 
range  accuracy  requirements.  Smaller  values  for  Performance  Level  (.75  to  .90) 
accommodate  the  realistic  possibility  that  under  some  circumstances  aircraft 
position  and  velocity  component  values  cannot  be  determined  accurately. 

Because  it  is  not  ideal  to  aggregate  the  accuracy  requirements  over  all 
scenarios  investigated,  the  individual  scenario  results  should  each  be  applied 
to  the  particular  portions  of  the  TACTS  range  on  vAiich  the  scenario  weapon  and 
target  are  to  be  exercised. 
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APPENDIX  A 


GENERAL  HETHCX)  FOR  DETERMINING 
TACTS  RANGE  ACCURACY  GUIDELINES 

The  general  method  for  determining  TACTS  range  accuracy  guidelines  is 
presented  in  this  appendix.  The  method  inplements  the  following  goal: 

Estimate  the  amounts  by  vAiich  aircraft  position  and  velocity 
ccmponent  values,  as  determined  by  TACTS  range  tel^netry  and 
processing,  could  be  in  error  while  requiring  that  TACTS  weapon 
simulations  compute  a  reasonably  accurate  value  for  Measure  of 
Effectiveness  (MOE),  reasonably  often. 

The  method  is  applicable  to  both  discrete  and  continuous  MOEs. 

Figure  A-1  depicts,  mathematically,  the  relationship  between  launch  param¬ 
eter  values  and  computed  MOE  (Score)  as  determined  by  a  TACTS  weapon 
simulation. 


A-1 


p  ;  L  -»  S  is  a  many-to-one  mapping  of  launch  conditions  into  scores.  This 
function  p  is  evalxxated  by  the  TACTS  weapon  simulations. 

Denote  /eh  and  5^ the  error  in  X  at  the  time  of  weapon  launch.  The  correct 
score  for  launch  condition  i  is  given  by  p(JI).  The  computed  score  resulting 
from  the  TACTS  range-measured  data  Jl  +  SX  is  p(/  +  5^).  Conditions  that  these 
two  scores  are  close  may  be  written 

|p(<+5^)  -  pU)|<e 
or 

h(jl+SJ)  -  p(^)|<e  •  p(i) 

in  which  eOO)  is  a  prescribed  tolerance.  The  first  of  the  conditions  places 
ein  eibsolute  constraint  on  accuracy,  independent  of  X  ;  the  second  specifies  a 
relative  accuracy  (i.e.,  p(  +  6J()  is  to  differ  from  p(j()  by  no  more  than  lOOe 
percent).  Note  that  if  the  scores  p(  Sf)  and  p(Jl)  must  match  exactly,  then 
choose  e  «  0.  Condition  (A-1)  is  probably  more  appropriate  for  discrete  noes 
(such  as  crater  diameter,  number  of  hits,  etc.);  the  second  condition  (A-2)  is 
more  appropriate  for  continuous  MOEs  (such  as  Pj^,  miss  distance,  percent  of 
target  destroyed,  etc.). 

Denote  the  absolute  difference  between  p( /  +  8^)  and  p(i)  by 

D(j(,5^)  -  |p(i  +  5/)  -  p(^)|  (A-3) 

Then  the  statistical  requirement  for  TACTS  range  accuracy  may  be  written 

Pr  {D(;?,S/)^e}  »  .95  (A-4) 

or 

Pr  (D(l6j)<e^p(j)}  >  .95  (A-5) 

Note  that  the  .95  is  chosen  arbitrarily  for  illustration. 


(A-1) 

(A-2) 


A-2 


Let  be  the  joint  probability  density  function  for  launch  conditions^ 

euid  the  TACTS  reuage  error  in  lamch  conditions  5^.  Requirements  (A-4,  -5)  are 
equivalent  to 

Jj^f(i,SjJ)did(5i)  ^  .95  {A-6) 

in  >diich  region  R  »  ^  e} 

or  alternately,  R  »  0  ^  e  p(J?)}. 

The  value  of  the  integral  in  {A-6)  depends  upon  f  in  particular,  if 

(Tjl  »  0  (and  -  0  so  that  TACTS  range  errors  are  unbiased),  then  the  value 
of  the  integral  is  1;  therefore,  all  scores  are  correctly  determined  by  the 
weapon  simulation  because  all  input  is  correct. 


As  increases,  the  integral  in  (A-6)  decreases  in  value.  There  is  a 
maximum  value  .of  <Tjjj  for  vdiich  condition  (A— 6)  holds.  For  the  set  of  launch 
platfonn/pilot/target/environmental  specifics  represented  in  f (.,.),  that  maxi¬ 
mum  (jj  is  the  TACTS  range  accuracy  required  to  guarantee  condition  (A-6). 

C 


As  a  practical  matter  f(.,.)  cannot  be  known,  so  a  reasonable  approximation 
must  be  employed  for  the  determination  of  the  integral  in  (A— 6)  as  a  function 
of  Or  more  correctly,  a  reasonable  approximation  for  the  integral  in 
(A-6)  as  a  function  of  SA  niust  be  used,  and  from  that  approximation,  the  value 
of  determined. 

One  such  method  of  approximation  is  as  follows;  let  jfeL  be  a  typical  set 
of  actual  launch  conditions.  Such  actual  launch  conditions  represent  a 
realistic  variation  on  intended  launch  conditions  recommended  by  tactical 
manuals  and  other  sources  of  tactical  guidance.  Evaluate  p(J).  Then  in  one 
conponent  of  >(  at  a  time  find  the  largest  and  smallest  variations  8J  for  which 
D(A»5j|)  ^  s  tor  D(j(,5i)  <  e  p(J{)]. 


For  several  values  of  i  representing  realistic  variations  in  launch  con¬ 
ditions,  determine  lower  and  upper  bounds  on  each  of  the  position  and  velocity 
parameter  values.  The  lower  bound  is  the  amount  by  which  the  parameter  can  be 
decreased  without  changing  the  value  of  p  substantially  -  similarly  for  the 
upper  bound. 

For  convenience  denote  the  parameters  by  i  »  1,...6,  with  «  x,  >^2  “ 
y.  ^3  -  ^4  -  ^5  -  Vy,  ig  -  v^.  Let  the  lower  and  upper  bounds  he  if 

and  Let  values  for  the  r®  set  of  launch  conditions  be  denoted  i^~(r)  , 

i.-"(rh  i  -  1,...,6. 

Suppose  that  the  TACTS  range  standard  deviations  for  the  determination  of 
these  parameter  values  are  i  «  1, . . . .6. 

If  the  errors  are  unbiased,  normal,  and  independent,  then 

M(r)  »  Pr  {D(j(,Si)  <  e}  (or  Pr{D(j?,5jJ)  ^  p(.^)}) 


(A-7) 


in  which  0  is  the  cumulative  normal  distribution  function. 


The  factor 


I 


is  the  probability  that  the  TACTS 


range  determines  a  launch  platform  x-coordinate  value  between  (r)  andy^j^^(r) 
and  similarly  for  other  position  and  velocity  conqponent  values.  The  product 
determines  the  joint  probcdsility  that  all  6  values  are  determined  with 
sufficient  accuracy. 


A-4 


If  the  conpated  simulation  score  is  to  be  accurate  95  percent  of  the  time, 
then  set 


Nr 

PL  =  - y  M(r)  -  .95  (A-8) 

Nr  ^ 
r-1 

in  which  Nr  -  number  of  simulation  runs,  PL  »  Performance  Level  ( .95  is  for 
illustration  only),  and  M(r)  is  given  in  (A-7).  The  function  PL  decreases  with 
increasing  Search  (numerically)  for  combinations  of  for  which  (A-8) 

holds.  These  combinations  of  values  specify  accuracy  required  on  the  TACTS 
range.  J^pendix  C  presents  the  listings  of  BASIC  programs  with  the  numerical 
procedures  for  conpating  desired  combinations  of  <t^'s  from  weapon  simulation 
data. 


A-5 


APPENDIX  B 

ESTIMATION  OF  VARIANCES  IN  LAUNCH  PARAMETERS 
FROM  MEASURES  OF  ACCURACY 


In  order  to  simulate  actual  launch  parameter  values,  estimates  for  pilot  ability  to  attain  in¬ 
tended  launch  parameters  must  be  available.  This  Appendix  describes  a  method  of  estimation 
employed  in  the  TACTS  weapon  simulations  for  purposes  of  this  analysis.  The  method  is  described 
in  the  context  of  bombing  a  fixed  target,  and  later  applied  in  air-to-air  gunnery. 

Denote: 

W.=  {x,y,z,9,tj),s) 

2C=(x,3/,z,I4,V„V;), 

each  a  vector  of  launch  parameters  at  the  time  of  bomb  release.  Coordinate  x  is  the  horizontal 
distance  between  the  aircraft  and  target;  y  is  the  aircraft’s  horizontal  position  perpendicular  to  the 
X-axis;  and  z  is  .the  vertical  distance  between  the  aircraft  and  the  target.  Dynamic  parameters  9,  rp, 
and  s  are  aircraft  pitch,  heading,  and  speed  respectively,  and  velocity  components  V^,  Vy,  and  V* 
are  determined  by  values  for  9,  and  s,  and  conversely.  A  pilot’s  ability  to  attain  intended  launch 
parameter  values  is  represented  by  the  variances  in  the  differences  between  actual  and  intended 
values. 

Let  AXh,  and  AY/i  be  the  differences  between  actual  and  intended  hit  position  coordinates. 
If  intended  launch  parameter  values  are  attained,  then  AXh  and  AYh  will  each  be  zero.  But  in 
general. 


,  dXn^  ,  dXh^  ,  dXn^  .  dXn^,,  ,  dX^ 


,  +  ^AV, 


dv. 


AYh  «  -I- 

OX 


dYh.  ^dYh.  ^dYh..,^dYh.,, 

■Ay  -I-  -^Az  +  -^AV;  +  TTrAVy 


dy 


dz 


dVr 


dv. 


•t".. 


in  which  AX.  =  (Ax,  Ay,  Az,  AVx,  AVy,  AVz)  are  the  errors  in  X_  at  launch. 

The  method  which  follows  estimates  variances  for  the  errors  in  AX.  subject  to  the  “known” 
conditions  on  the  hit  position  errors  AXh,  AYh.  The  “known”  conditions  are  that  bombing  accuracy 
is  “5  mils”  and  gunning  accuracy  is  “14  mils”. 

As  a  first  step  the  partial  derivatives  . . . ,  axe  estimated.  Let  Wj  denote  a  set  of 
intended  launch  parameter  values  as  determined  by  tactical  guidance.  Run  the  bomb  simulation  to 
obtain  bomb  impact  position  -  call  it  (Xh,Yh)  -  to  be  thought  of  as  the  intended  hit  position. 

Then  for  each  component  in  2C  (one  at  a  time),  increase  the  parameter  by  a  small  amount  and 
rerun  the  bomb  simulation  with  the  perturbed  launch  parameter  value.  (Perturbations  in  Vx,  Vy, 
Vz  must  be  converted  to  perturbations  in  9,  tp,  and  s  before  running).  Obtain  the  bomb  impact 
positions  under  these  six  variations.  (Variations  in  x  and  y  need  not  be  run  because  the  bomb 
trajectory  is  simply  translated  so  that  miss  distances  on  the  ground  are  exactly  the  values  in  x  and 
y  variation.) 

Denote  the  changes  in  impact  positions  {Xh,Yh)  resulting  from  parameter  variations  by 


AXh{Ax)[=  Ax],  Ay/,(Ax)[=  0],  AXA(Ay)[=  0],AYh(Ay)[=  Ay], 

AXhiAz),AYh{Az)[=  0],  AXh(AV^),  An(AF,)[=  0],  (B-1) 

AX/.(AF,)[«  0],An(AV,),  AXh{AVz),AYhiAVz)[=  0]. 

Values  known  in  advance  are  given  in  “[  ]”. 
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Then  as  a  first-order  approximation, 


An«A,+  {^^^}.AV. 


(B-2) 


in  which  the  terms  in  “{  }”  are  the  numerical  estimates  of  partial  derivatives.  Alternative  numer¬ 

ical  estimation  formulas  for  partial  derivatives  may  be  substituted  for  those  used  to  obtain  more 
accurate  estimates. 

To  simplify  notation,  abbreviate 


and 


rAr,(Av„)i 
1  Av;  /  - 


.If  it  is  assumed  that  the  errors  A2C  =  (Aar,  Ay,  Az,  A V.,  A V„,  AV,)  in  X  at  launch  are  independent 
normally  distributed  random  variables  with  0  means,  then  =  J?(Ax*],  =  £^[Ay*],  = 

r[AiAy](=  0),  etc.,  in  which  i?[-)  is  the  “expectation”  function.  Consequently,  Axl,  aV/  in 
equation  (B-2)  are  normal  and  independent. 


e[axh]  =  i;[An]  =  0, 


and 

<^AX^  «  4.  +  •  <Tiv.  +  DV,^  .  (riv; 

•  4 v;  • 

Suppose  that  bombing  accuracy  is  5  mils.  Then  the  median  miss  distance  in  the  plane  perpen¬ 
dicular  to  the  weapon  trajectory  at  the  point  of  impact  is  5L  feet  as  shown  in  the  following  Figure, 
where  L  =  trajectory  length  in  thousands  of  feet  (e.g,,  I  =*  2  for  a  2000  ft  trajectory). 


z 


The  elliptical  area  on  the  ground  plane  which  contains  50  percent  of  the  bomb  hits  has  axes 


Ay  =  5L  and  Ag 


51 

sin  9/ 


in  which  0/  =  weapon  impact  angle.  Thus 


~  ~  sin©/' 
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(B-4) 


The  probability  within  a  circle  of  radius  r  centered  at  the  mean  for  a  circular  normal  random 
variable  (17,  V)  with  <tu  =  cry  =  o-  is  given  by  1  —  exp(—r^/2<T^).  If  this  has  value  .5,  then  the 
corresponding  r  =  <t\/2  •  In  2  is  the  median  miss  distance,  where  v^2  •  ln2  «  1.177. 

Thus  for  CEP  =  5X,  the  corresponding  standard  deviation  in  the  plane  perpendicular  to  the 
trajectory  at  impact  is  <r  =  5X/1.177.  Thus 


= 


5i 

1.177 


and  = 


1 


5X 


1.177  sin0/’ 


Now  abbreviate: 


C.= 


5i 


1.177  sin0/ 

C  =-^ 

1.177' 

Then  equations  (B-3)  and  (B-5)  together  imply  the  constraints 

4,  +  Dz^  .  4,  +  DV^  .  aly^  +  DV^  •  cj 


(B-5) 


+  DV^-<’iv.=C, 


(B-6) 


on  the  six  error  standard  deviations. 

Under  the  normality  and  independence  assumptions  on  the  components  if  AX.,  the  joint  prob¬ 
ability  density  function  for  the  six  error  components  is  proportional  to  Q: 


Q  =  exp 


<^AV, 


)  +( 


(B-7) 


Among  all  volumes  which  contain  the  six  error  values  (Aa:,Ay,  Az,  AF®,  AFy,  AF^)  50  percent  of 
the  time  (so  that  the  resulting  bombing  accuracy  is  5  mils),  the  smallest  volume  is  defined  by 
Q  >  Ki  for  some  constant  K\  (whose  actual  value  is  irrelevant). 

To  confirm  that  Q  >  K\  defines  the  smallest  volume,  note  that  f(x)(=  /(xi, . .  .,x„))  is  a 
probability  density  function  for  random  variable  X  if  and  only  if  /(x)  >  0  for  all  £  and  f(x)dV  = 
1.  E’^  is  all  of  n-dimensional  space  (=  {(xi, . .  .,x„)  |  — oo  <  Xi  <  -foo,  i  =  1, . .  .,n})  and  dV  is  the 
volume  differential  (=  dx\dx2  •  •  •  dx„). 

A  subset  B  C  IS”  is  a  confidence  volume  for  X  at  confidence  level  C,  0  <  C  <  1,  if  Jg  f{x)dV  = 
C.  Define 

A{k)  =  {x  |/(x)  >  fc). 

Then  A(0)  =  E'^  and  A(k  +  e)  C  A{k)  for  all  c  >  0.  As  k  increases,  A{k)  gets  smaller;that  is 
limfc_co  A(A;)  =  0,  the  empty  set:  further, 

Pr(A(^'))  =  /  fix)dV  >  f  kdV  =  k  [  dV  =  k’  Vol{Aik)) 

JA{k)  JAik)  JAik) 

and  because  Pr(A(A:))  <  1,  then  limfc_»oo  k  •  Vol(A{k))  <  1  and  therefore  limfc_oo  Vol(A(k))  =  0. 
Thus  limfc-^oo  Pr(-^(^))  =  0.  Consequently,  for  fixed  C,  0  <  C  <  1,  there  is  a  constant  k  =  k{C) 
such  that  provided  that  /  is  continuous  with  non-zero  partial  derivatives. 

The  set  A(fc(C'))  (abbreviated  A)  has  smallest  volume  of  all  confidence  volumes  at  level  C,  as 
shown  by  the  following  (B  is  any  other  confidence  volume  at  level  C): 

/  /(x)dF  =  /  f{x)dV  +  [  f(x)dV  =  C, 

Ja  JAnB  JAnS 
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so  that  by  subtraction, 


so  adding, 


/  fix)dV  =  /  f(x)dV  +  /  f{x)dV  =  C, 

Jb  Jahb  JAnB 

/  fix)dV=  f  f(x)dV. 

JAnB  JAnB 

f  f{x)dv>  [  kdV  =  kVol(AnB) 

JAnB  JAnB 

f  f{x)dV<  f  kdV  =  kVol(AnB). 

JAnB  JAnB 

Vol{A  nB)<  Vol(A  n  B), 
Vol{AnB)  =  Vol{Af]B), 

Vol{A)  <  Vol(B). 


The  conclusion  is  that  of  all  confidence  volumes  at  confidence  level  C,  the  one  of  smallest  volume 
is  A(k(C)). 

The  shape  of  that  smallest  volume  is  an  ellipsoid.  Set  Q  >  Ki  {Q  is  defined  in  equation  (B-7)) 
and  take  the  natural  logarithm,  obtaining  the  volume  defined  by 

+  {^)  +  +  {^)  <  K] 

(TA®  <7Ay  (f&.z  (TAV,  O^AVi  o'AV, 

for  some  constant  This  volume  is  an  ellipsoid  in  the  6-dimensional  space  of  errors  (whose  points 
are  labeled  (Ax,  Ay,  Az,  AVa:>  AFy,  AV^)),  with  semi-axes  proportional  to  o"Ay>  «7-az, 

(TaKv.  and  ctav, • 

The  volume  of  A{k{C))  is  proportional  to  the  product  of  the  ellipsoid  semi-axes,  namely; 


V  —  K3  ’  <TAx<XAy<JAz<XAV^(JAVy(JAV, 


(B-8) 


for  a  constant  K3  (whose  value  is  irrelevant). 

The  two  constraints  of  equation  (B-6)  do  not  provide  sufficient  information  for  the  determi¬ 
nation  of  the  six  error  standard  deviations.  Further,  specific  information  regarding  relative  sizes 
of  the  six  error  standard  deviations  is  not  available.  Consequently,  a  mathematical  optimization  is 
employed  in  which  values  for  (o'Aa:jO'Ay»o’Az»o‘AV,  jO’av^jO'av,  )  are  chosen  to  maximize  the  volume 
in  equation  (B-8)  subject  to  the  constraints  of  equation  (B-6).  The  optimization  is  performed 
analytically  using  the  method  of  Lagrange  multipliers.  Notice  that  the  optimization 


maJtF  =  m^/V3  ’  (JAx<JAy<XAz(XAV^<JAVyCrAV, 

{<T'4}  {<7'3} 


is  equivalent  to  the  two  optimizations 


max  <Jax<Jaz(Xav^<Jav, 

{O’ Ax  .«^A«  .‘^AV’i  I'^AV,  } 

S.T.  4,  -1-  L>z2  .  al,  -h  DV^  •  4v;  + 
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and 


max 

{<r^y,<7AVy  } 

S.T.  (Tly  +  DV;-<Tly^=Cl 

because  the  objective  function  V  and  the  constraints  separate  into  two  independent  sets. 

Each  of  these  optimizations  has  the  form 

n 

maxTT  Xi 

S.T. 

t=l 

which  is  a  constrained  optimization  problem.  This  can  be  solved  by  maximizing  the  Lagrangian 
function 


Differentiation  yields 


i=l  t=l 

|^=  n  -  2An?X.- =  0,  f=l,...,n 


j=i 

dX 


*=1 


The  first  yields  (multiplying  by  Xi) 


]JXj  =  2XaJX?,  i  =  l,...,n 

i=i 

so  that  aU  terms  a^Xf  are  equal.  Let  their  common  value  be  u^;  then  from  the  condition  —  =  0 
above, 

—  I?  =  0,  nv?  =  b^,  u  = 

h.  ’  ’ 


and  consequently 


^i  = 


\b\ 


\ai\y/n 


^  ^  n. 


Apply  this  general  result  to  the  first  optimization  above  by  letting  Xi  =  ^2  ==  <^A2j 

-^3  =  <^AVr  >  ^4  =  «^AV, ,  of  =  1,  =  Dz^,  03  =  DVx^ ,  =  DVj^ ,  6^  =  cl,  and  n  =  4.  Similarly, 

apply  this  general  result  to  the  second  optimization  above  by  letting  Xi  =  cr^y,  X2  =  <tav,,  a?  =  1» 
02  =  DVy^,  =  Cl ,  and  n  =  2.  As  a  result,  the  following  estimates  are  obtained: 


Cx  ^  Cy  ^  _  Cx 

^^Ax-Y^  '^^^-2\Dzy 


=  TTTTTfTr;  <^AV,  = 


2\DVxy 


V2\DVy\' 


<TaV,  = 


2|T»v; 


(B-9a) 

(B-9b) 
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These  error  standard  deviation  estimates  in  equations  (B-9)  are  used  to  simulate  actual  launch 
conditions  with  respect  to  stationary  targets. 

A  similar  approach  is  used  to  determine  variances  in  launch  parameter  values  for  air-to-air 
gunnery. 

Assume  that  the  Closest  Point  of  Approach  (CPA)  occurs  in  the  r-s  plane  perpendicular  to 
Uic  bullet  s  path.  At  CPA  the  target  is  at  position  (0,r,s)  as  shown  in  the  following  figure,  where 
CPA  =  £,  and  -f  . 


target  s 


Assume  at  CPA  that  r,  s  are  values  of  independent,  normal  random  variables,  each  with  means 
M  =  0  and  equal  variances  cr^.  Then  (£/<r)2  is  a  random  variable  with  Chi-square  distribution  with 
2  degrees  of  freedom  (abbreviated  2d/).  The  density  function  for  tliis  random  variable  is 

1  X 

If  m  is  the  median  miss  distance,  then  Pr(£  <  m)  =  .5;  that  is,  Pr(£2  <  m^)  =  .5  or  < 

(m/cr)  )  =  .5.  Thus  {m/cry  is  the  median  M  for  the  random  variable  with  2df.  The  numerical 
value  AI  is  determined  as  follows;  if 


then 


I 


M. 


M 


|exp(-|)da:  =  .5, 


1  -  cxp(- y)  =  .5  and  Af  =  ln4  «  1.386. 


The  Air  Comba,t/MissiIe  EvaluAtion  Study,  from  1976  (reference  c)  suggests  that  m  =  h\L  fa 
“14-mil  error”)  so  that 


(Mi)  =,„4  = 

\  J  ln4  ’ 


where  is  the  common  variance  for  r  and  s. 

Next,  determine  the  variance  for  £; 

From  E[{^f]  =  2  (the  expected  value  of  a  random  variable  =  df)  there  follows 

=  2(T^  (using  linearity  of  the  expectation  function). 


(B-10) 

(B-11) 

(B-12) 
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Further,  =  E[l].  Let  w  =.(£/<r)^;  then  tz>  is  a  random  variable  with  2d/;  t  =  <Ty/w;  and 
=  /o°°  cry/wfy2{w)dw  which  can  be  simplified  to  yield 


=  |v^* 

Putting  results  in  equations  (B-10)  through  (B-13)  together  yields 


(B-13) 


a  I  =  2a^  — 


2^(7^  4  —  X 


and 


^2  ..^2  4-7r(14X)2 

2  ln4  ’ 


(B-14) 


in  which  C}  is  an  abbreviation  for  the  constant  given.  Standard  deviation  at  in  niiss  distance  I  is 
a  function  of  deviations  in  parameters  x,y,z,Vx,Vy^Vx  at  the  time  of  gun  firing;  the  actual  miss 
distance  is  given  by 


A^ 


^Ax  +  -  Aj  +  gj  A.  +  —  Ar.  +  gj^AF,  +  g-An. 


Under  the  assumption  that  the  above  parameter  deviations  AX  are  unbiased  and  independent, 


a\  «  E[^l\ 


As  in  the  treatment  of  stationary  targets,  abbreviate  the  estimate  for  ||  as  Dx,  . 
as  DVz.  Then  perform  the  following  optimization: 


(B-15) 
and  for 


max 

{a'j} 


<^AV„<^AV’, 


S.T.  Di"  •  „l,  +  •  4,  +  •  4,  +  DV!  •  4v.  +  DV^  •  4^.,  +  DVl  ■  4^-,  =  C?, 


obtaining  the  estimates 


S  |g'l  s  ICtl  ~  |C<| 

"  |Di|V6’  “  |Dx|v^’ 


(B-16a) 


S'av; 


I-DT4IV6’  |X>V;|V6’  |X>U^|a/6' 


(B-16b) 


If  an  estimate  |Z?p|  for  a  partial  derivative  |^|  for  a  parameter  p  is  nearly  0  (as  is  the  case 
for  parameters  x  and  Vx),  then  t  is  insensitive  to  p.  In  this  instance,  we  may  take  dp  =  0,  because 
the  dp's  are  employed  in  simulating  errors  in  p  at  the  time  of  gun  firing,  and  the  error  in  p  would 
be  inconsequential  insofar  as  simulating  bullet  CPA  and  measures  of  effectiveness  depending  upon 
CPA. 

These  error  standard  deviation  estimates  in  equations  (B-16)  are  used  to  simulate  actual  launch 
conditions  with  respect  to  airborne  targets. 
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APPENDIX  C 


PROGRAM  LISTINGS 


Presented  in  this  appendix  are  listings  of  two  BASIC  programs  used  in  the 
determination  of  TACTS  range  accuracy  requirements. 

Program  DISCRETE  FIT  implements  the  methodology  of  Appendix  A  for  the 
computations  of  position  and  velocity  component  standard  deviations.  Equations 
{A-7)  and  (A-8)  are  involved,  together  with  data  produced  by  TACTS  weapon 
simulation  runs. 

The  other  program,  DFIT_BAT_PLOT  performs  the  computations  of  DISCRETE_FIT 
automatically  for  a  variety  of  Performance  Levels  and  data  frcxn  one  or  several 
TACTS  weapon  simulation  nans.  Analyzed  results  are  stored  in  data  files  for 
subsequent  graphics  presentations. 
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DISCRETE  -  FIT 


29  Jan  1987  10:31:09 

100  REn  DISCRETE_FIT 
110  ! 

120  REM  WRITTEN:  AUGUST  14,  1986 

130  REM  LAST  MODIFIED:  SEPTEMBER  4,  1906 
140  .  ! 

150  COM  /Consta/  Pl.BiS) 

160  COM  /Ratios/  Ratio(l:6) 

170  COM  /Bounds/  Lower, Upper 
180  COM  /Output/  Sig(l:6) 

190  COM  /Intervals/  Nruns ,X( 1 ,6 . 100 ) 

200  DIM  Si9i(l:6) ,Nin$(801 
210  ! 

220  PRINT  PAGE 

230  GOSUB  350  !  INITIALIZATION 

240  GOSUB  660  !  READ  SIMULATION  DATA 

250  GOSUB  860  !  INPUT  PRECfSION 

260  GOSUB  950  .!  PRINT  HEADER 

270  WHILE  FNPleve I  CO ,Siqv) >Desp I 

280  GOSUB  1140!  BISECTION  SEARCH 

290  PRINT  USING  "  6  C  2X  ,  4D .  2D  ,X )  *' jS  ig  C  *  ) 

300  Sigv“Sigv+Oe Is igv 
310  END  WHILE 
320  PRINTER  IS  CRT 
330  END 
340  ! 

350  REM  initialization 
360  ! 

370  DATA  .2316419 

380  data  0, .31938153 .-.356563782,1.781477937,-1.821255978 ,1.330274429 

390  READ  PI 

400  FOR  Param*0  TO  5 

410  READ  BCParam) 

420  NEXT  Param 
430  I 

440  Lower>0 
450  Upper*! 

460  ! 

470  Eps t Ion* . 00001 
480  Sigv»0 
490  Delsigv*.01 
500  ! 

510  SigSCl)*"  SIGMA-X 
520  Sig«(2)«"  SIGMA-Y  " 

530  Sig$(3)-"  SIGMA-Z  “ 

540  Sig«(4)-"  SIGMA-Ux  " 

550  Sig*(5)-"  SIGMA-Uy  " 

560  Sig«C6)-"  SIGMA-Uz  " 

570  ! 

500  RatioCD-l 
590  Ratio(2)-l 
600  Ratio(3)*4 
610  Ratio(4)«l 
620  Rat  10  C5 )■! 

630  Ratio(6)*4 
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640  RETURN 


650 
660 
670 
680 
690 
700 
710 
720 
730 
740 
750 
760 
770 
780 
790 
800 
810 
820 
830 
840 
850 
360 
870 
880 
890 
900 
910 
920 
930 
940 
950 
960 
970 
980 
990 
1000 
1010 
1020 
1030 
1040 
1050 
1060 
1070 
1080 
1090 
1100 
1110 
112  0 
1130 
1140 
1150 
1160 
1170 
1180 
1190 


REM  READ  SIMULATION  DATA 

j 

REM  X(l,6,100)  I  1  •  Lower, Upper 

!  6  •  X,Y,Z,Ux ,Uy ,Uz 
!  100  •  Runs 


PRINT  "  INPUT  NAME  OF  DATAFILE? 
INPUT  "’^Ninl 
PRINT  Nin$ 

PRINT 

ASSIGN  *1  TO  Nin» 

READ  ♦IjNruns 
FOR  R-1  TO  Nruns 
FOR  Pa ram- 1  TO  6 

READ  ♦!  JX(  0  ,Param,R)  ,X(  1  ,Paratn  ,R) 
NEXT  Pa ram 
NEXT  R 

ASSIGN  ♦!  TO  * 

RETURN 


REM  INPUT  PRECISION 


PRINT 

PRINT  "INPUT  DESIRED  PERFORMANCE  LEUEL  (Between  0  &  1 )  -  J 
INPUT  "",DeapI 
PRINT  Dcspl 
PRINT 
RETURN 
! 

REM  PRINT  HEADER 

» 

PRINT  "HARDCOPY  (Y  or  N): 

INPUT  "",Hc$ 

IP  Hc*»"Y"  THEN  PRINTER  IS  402 
PRINT  PAGE 

PRINT  “NAME  OF  DATAFILE?  ";Nin$ 

PR  1  NT 

PRINT  "DESIRED  PERFORMANCE  LEUEL  -  ";Oespl 

PRINT  .. 

PRINT  "  POSITION  ERROR  RATIO  -  X?Y:2  -  » 

PRINT  USING  "0,A,D,A,0"  ;Rat  lo  ( 1)  I  •' ?  "  }Rat  io  (2 ) } ’[  s  “  ;Ratio(3. 
PRINT  "  UELOCITY  ERROR  RATIO  -  Ux?Uy? Uz  -  '  > 

PRINT  USING  "0,A,D,A,D" jRat lo ( 4 ) } " : " ;Rat lo ( 5  ) ; " ? " }Rat lo (6 . 

PRINT 

PRINT  USING  "2X,6(10A)" ;Sig»(* ) 

PRINT 

RETURN 

j 

rIm  ^FNPLEUEL(^^*!'')*”'iS  a  DECREASING  FUNCTION  IN  BOTH  ARGUMENTS 


Sig 1 lo* . 01 

WHILE  FNPleve I (Sig 1 lo ,Sigv)<Desp I 
Sigl lo*Sigl Io/2 
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1200  END  WHILE 
1210  S»glhi»l 

1220  WHILE  FNP leva  I CSig Ih X ,Sigv) >Oe*p  I 
1230  Siglhi-2*Siglhi 
1240  END  WHILE 

1250  SigImid-.5*(Sigl lo>Siglhi ) 

1260  PImxd-FNPleve  l(Sigltnid,Sigv) 

1270  while  ABS(PImxd-Oesp U >Epa i Ion 
1200  IF  Plinid>Oespl  THEN 
1290  Sig I lo*Sig Imid 

1300  ELSE 

1310  Sig Ih i»Sxg Imid 

1320  END  IF 

1330  Sig Imid« .5*(Sig I lo^Sig Ihx ) 

1340  Plmid-FNPlevel (Siglmid,Sxgv) 

1350  END  WHILE 
1360  RETURN 
1370  ! 

1380  END 
1390  ! 

1400  DEF  FNPles^el  CSx  .Sv) 

1410  COM  /Bounds/  Lower , Upper 

1420  COM  /Ratios/  Ratio  I*) 

1430  con  /Output/  SigC*) 

1440  COM  /Intervals/  Nruns,X(*) 

1450  DIM  S(l:6) 

1460  Prob_sum“0 
1470  S(l)-Sx*Rat 10(1) 

1480  S(2)«Sx*Rat io(2 ) 

1490  S(3 )-Sx*Rat 10 (3 ) 

1500  S(4)»Sv*Rat 10(4) 

1510  S(5 )-Sv*Rat 10(5 ) 

1520  S(6 )«Sv*Rat 10(6 ) 

1530  FOR  Run«l  TO  Nruns 
1540  Prob"! 

1550  FOR  Param-l  TO  6 

1560  IF  S(Param)-0  THEN 

1570  Factor-1 

1580  ELSE  ^  ^ 

1590  Factor«FNP(X(Upper ,Param,Run )/S(Param) ) 

1600  Factor-Factor-FNP(X(Lower ,Param,Run)/S(Param) ) 

1610  END  IF 

1620  Prob-Prob*Factor 

1630  NEXT  Param 

1640  Prob_sum»Prob_sum*Prob 

1650  NEXT  Run 

1660  Pleve l*Prob_sum/Nruns 

1670  MAT  Siq-S 

1680  RETURN  P level 

1690  FNEND 

1700  ! 

1710  DEF  FN2CX) 

1720  IF  ABS(X)>5  THEN 

1730  Z«0 

1740  ELSE 

1750  2-(l/SQR(2*PI ) ) *EXP( -X*X/2 ) 
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1760  END  IF 
1770  RETURN  2 
1780  FNENO 
1790  ! 

1800  DEF  FNP(Xin) 

1810  COM  /Consts/  P1,BC*) 

1820  •  IF  Xin>0  THEN 

1830  X«Xin 

1840  ELSE 

1850  X»-Xin 

I860  END  IF 

1870  T-1/(1*P1*X) 

1880  Temp«0 

1890  FOR  I«5  TO  0  STEP  -1 
1900  Temp*T*Temp^B( I ) 

1910  NEXT  I 
1920  P«l-FN2CX)*Temp 
19"30  IF  Xin<0  THEN  P-l-P 

1940  RETURN  P 

1950  FNENO 


DFIT  -  BAT  -  PLOT 


29  Jan  1987 


10:33:02 


100 

110 

120 

130 

140 

150 

160 

170 

130 

190 

200 

210 

220 

230 

240 

250 

260 

270 

280 

290 

300 

310 

320 

330 

340 

350 

360 

370 

380 

390 

400 

410 

420 

430 

440 

450 

460 

470 

480 

490 

500 

510 

520 

530 

540 

550 

560 

570 

580 

590 

600 

610 

620 

630 


REM  0F1T_8AT_PL0T 

REM  URITTEN:  SEPTEMBER  11,  1986 

REM  LAST  MODIFIED:  OCTOBER  2.  1986 
! 

COM  /Consts/  P1,B(5) 

COM  /Ratios/  Ratto(l:6) 

COM  /Bound#/  Lower, Upper 
COM  /Output/  Stg(l:6) 

COM  /Interval#/  Nrun# ,Xt 1 ,6 . 100 ) 

DIM  Sig»ll:6), Index f i le$t30],Nin$C80] 

DIM  Xout(0;50),Yout(0:50) 

! 

PRINT  PAGE 

GOSUB  710  !  INITIALIZATION 

GOSUB  1010!  DETERMINE  CASES  TO  BE  RUN 

FOR  Ca#e»l  TO  Nca#e# 

GOSUB  1210!  READ  SIMULATION  DATA 
FOR  De#pl-.95  TO  .7  STEP  -.05 
GOSUB  1870!  BISECTION  SEARCH  FOR  Maxsigv 
Neva4-25  !  NUMBER  OF  EUALUATIONS  EQUALLY  SPACED 
Ext  pa-4  !  EXTRA  EUALUATIONS  TO  SMOOTH  UP  GRAPH 

De I# igv-Max# igv/Neva 1 
GOSUB  1390!  PRINT  HEADER 
FOR  Eva  1-0  TO  Neva  I ♦Extra 
SELECT  Eva  I 
CASE  -Neva  1 ♦Extra 
GOSUB  1790 
CASE  ELSE 
IF  Eva  I < Neva  I  THEN 

S i gv-Eva I *De 1 # i gv 

E>-SE  , 

S  igv-  (Neva  l-l)*Del#igv^t  Eva  1  ♦  1-Neva  1  )  *De  I#  igv'  «.Extra*l) 

END  IF 

GOSUB  1550!  BISECTION  SEARCH 
END  SELECT 

PRINT  USING  '’6t2X,4D.2D,X)"  ;Sig(*  J 
Xout (Eval )-Sig(4) 

Yout (Eva  I )*Sig ( 1 ) 

NEXT  Eval 

GOSUB  Putoutp lot  data 
IF  NOT  He  THEN  WAIT  1 
NEXT  Oe#pl 
NEXT  Case 
PRINT  PAGE 
PRINTER  IS  CRT 
END 

I 

Putoutp lot  data:  ! ESTABLISH  PLOTDATAF ILES 

Pname»-Ca#ename«(Ca#e)i"_"&UAL»(Rat  io (3  ))&••«" &UAL$ ( 100*Oe#p  1) 
CREATE  DATA  Pnamei,20 
ASSIGN  *1  TO  Pnamet 
PRINT  #1 ;Neval*Extra^l 
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640 
690 
660 
670 
680 
690 
700 
710 
720 
730 
740 
75  0 
760 
770 
780 
790 
800 
810 
820 
830 
840 
850 
860 
870 
880 
890 
900 
910 
920 
930 
940 
950 
960 
970 
980 
990 
1000 
1010 
1020 
1030 
1040 
1050 
1060 
1070 
1080 
1090 
1100 
1110 
1120 
1130 
1140 
1150 
1160 
1170 
1180 
1190 


FOR  Ev«l»0  TO  N«v«l»Extr8 
PRINT  ♦! ;Xout (Ev8l ) ,Yout CEv8l ) 

NEXT  Ev8 1 
ASSIGN  ♦!  TO  • 

j 

RETURN 

I 

REM  INITIALIZATION 

! 

DATA  i|”jt9jL5J,-.J5656J782,l.791477'937,-1.82125»978,l.JJ027A429 

READ  PI 

FOR  Par«m-0  TO  5 
READ  8 (Pa ram) 

NEXT  Pa ram 

! 

Louar«0 

Upper«l 

I 

Eps  1  Ion* . 00001 
Del»igv*.05 

I 

Sig»(l)-"  SIGMA-X 
Sig*(2)-"  SIGMA-Y  " 

Sig|(3)-"  SIGMA-Z  " 

Sigt(4)-"  SIGMA-Ux  ■' 

Sig»(5)-"  SIGMA-Uy  ” 

Sig»(6)-"  SIGMA-Oz  ” 

I 

Ratio(l)-l 
Rat  io(2)-l 
Ratio(3)*l 
Ratio(4)*l 
Ra  t i o ( 5 ) ■  1 
Ratio(6)*l 
RETURN 
! 

REM  DETERNINE  CASES  TO  BE  RUN 

PRINT  "INPUT  INDEX  FILE  NAME  OF  CASES  TO  BE  RUN: 

INPUT  ,  Indaxf  i  la* 

PRINT  Indaxfilaf 
ASSIGN  ♦!  TO  Indaxfila* 

READ  tllNcasas 
FOR  Caaa-1  TO  Ncaaaa 
READ  #l;Caaanamat(Caaa) 

NEXT  Casa 
ASSIGN  #1  TO  * 

I 

PRINT  "HARDCOPY  (Y  or  N):  "; 

INPUT 
PRINT  He* 

Hc-(Hc»-"Y") 

IF  He  THEN  PRINTER  IS  401 

! 

RETURN 
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1200 

1210 

1220 

1230 

1240 

1250 

1260 

1270 

1280 

1290 

1300 

1310 

1320 

1330 

1340 

1350 

1360 

1370 

1380 

1390 

1400 

1410 

1420 

1430 

1440 

145  0 

146  0 
1470 
1480 
149  0 
1500 
1510 
1520 
1530 
1540 
1550 
1560 
1570 
1580 
1590 
1600 
1610 
1620 
1630 
1640 
1650 
1660 
1670 
1680 
1690 
1700 
1710 
1720 
173  0 
1740 
1750 


REM  READ  SIMULATION  DATA 

I 

REM  X(l,6,100)  !  1  •  Lower ,Upp«r 

!  6  -  X,Y,2,Ux,Uy,Uz 
I  100  ■  Runs 

! 

Nin*«C«senamet (Case ) 

PRINT 

ASSIGN  «1  TO  Nin$ 

READ  tUNruns 
FOR  R-l  TO  Nruns 
FOR  Parem-l  TO  6 

READ  ♦! ;X(0 ,Param,R) ,X( 1 ,Param,R) 
NEXT  Pa ram 
NEXT  R 

ASSIGN  #1  TO  * 

RETURN 


REM  PRINT  HEADER 


PR  I  NT  PAGE 

PRINT  "NAME  OF  DATAFILE:  " ;N»nt 
PR  I  NT 

PRINT  "DESIRED  PERFORMANCE  LEUEL  •  ";Oespl 

print  ^ 

PRINT  "  POSITION  ERROR  RATIO  -  X:Y:2  - 

PRINT  USING  "D ,  A,0,  A,0*’  ;Rat  xo  ( 1 )  J "  !  "  jRat  io  (2  )  ;  '•  s  "  ;Rat  lo  (3  ) 

PRINT  "  UELOCITY  ERROR  RATIO  -  UxsUyrUz  -  "J 

PRINT  USING  "D,A,D,A,D" ;Rat io(4) " ;Rat io(5 ) j 
PRINT 

PRINT  USING  "2X,6(10A)" ;Sigf(*) 

PRINT 

RETURN 

! 

REM  ^FNpSiEL(?^?*)^*i's  A  DECREASING  FUNCTION  IN  BOTH  ARGUMENTS 


;Rat 10 (6 ) 


Sigl  lo-. 01 

WHILE  FNPIeve 1 (Sigl  lo ,Sigv)<Oesp I 
Sigl lo-Sigllo/? 

END  WHILE 


Siglhi-1 

WHILE  FNPIevel (Siglhx ,Sigv)>Despl 
Siglhi-2«Siglhi 
END  WHILE 

Siglmid-.5«(SxgI Io*Siglhi ) 
PImid-FNPIevel (SiglmidjSxgv) 
WHILE  ABS(Plmid-Despl)>Epaxlon 
IF  Plinxd>Oeapl  THEN 
Sigl  lo«Sigltnxd 
ELSE 

Sxglhx*Siglmid 
END  IF 

Siglmid-.5»(Sigl  lo>*-SigIhi ) 
Plmxd"FNPleve l(5’almid,Sxgv) 
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1760 
1770 
1780 
1790 
1800 
1810 
1820 
1830 
1840 
1850 
1860 
1870 
1880 
1890 
1900 
1910 
1920 
1930 
1940 
1950 
1960 
1970 
1980 
1990 
2000 
2010 
2020 
2030 
2040 
2050 
2060 
2070 
2080 
2090 
2100 
2110 
2120 
2130 
2140 
2150 
2160 
2170 
2180 
2190 
2200 
2210 
2220 
223  0 
2240 
2250 
7260 
2270 
2280 
229  0 
2300 
2310 


END  UHILE 
RETURN 

! 

REM  SIS(NEUAL)  UALUES 

! 

MAT  Sig-(0> 

FOR  1-4  TO  6 
Sig(I)-Maxs  igv-Rat  lO  (.  I) 
NEXT  I 
RETURN 

I 

REM  BISECTION  SEAI 
REM  FNPLEUEL(.,. 


Sigl lo- . 01 
UHILE  FNPleveKO,! 

Sigl lo-Sig I lo/2 
END  UHILE 
Siglhi-1 

UHILE  FNPlevel(0,S 
Siqlhi-2*Siglhi 
END* UHILE 
Siglmid-.5*(Sigllc 
Plmid-FNPlevel  (0,£ 
UHILE  ABS(Plmid-Oa 
IF  P-lmid>Oesp  1  Tl- 
Sigl lo-Sig Imid 
ELSE 

Siglhi-Siglmxd 
END  IF 

Siglmxd-.5*(SigI I 
PImid-FNPleval (0 , 

END  UHILE 

Maxa igv-Sig Imid-Ep 

RETURN 

END 


ARGUMENTS 


naf  A  nue 


N.A.D.C. 


U  890  0069 


Tech.  Mem.  88-202-5 


DATE  DUE 


BORROWER’S  NAME 


OEF  FNPleveKSxjS*. 

COM  /Bounds/  Lows 
COM  /Ratios/  Rst i 
COM  /Output/  Sig( 

COM  /Intervals/  N 
DIM  S(l:6) 

Ppob_sum-0 
S(1 )-Sx*Rat io(l) 
S(2)-Sx*Rat io(2) 
S(3)-Sx*Ratio(3) 
S(4)-Sv*Rat io(4) 
S(5)-Sv*Rat ioC5) 
S(6)-Sv»Rat 10(6) 

FOR  Run-1  TO  Nruns 
ppob-1 

FOR  Param-l  TO  6 
IF  S(Papam)-0  THEN 
Factor-1 
ELSE 
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2320  F«ctor«FNP(X(Upp«r ,P«ram,Run )/S(Param; ) 

2330  Factor«Factor-FNP(X(Low8r  ,Param,Run)/'S(Param) ) 

2340  END  IF 

2350  Prob«Pr'ob*Factor 

2360  NEXT  Param 

2370  Prob_sum»PrQb_sum-»Prob 

2380  NEXT  Run 

2390  Pleva  l*Prob_sum/'Nruns 

2400  MAT  Sig-S 

2410  RETURN  Plaval 

2420  FNENO 

2430  ! 

2440  DEF  FN2(X) 

2450  IF  AeS(X)>5  THEN 

2460  2-0  r  ' 

2470  ELSE  i  r  ^ 

2480  2-C1/SQR(2*PI) )*EXPC-X*X/2) 

2490  END  IF  " 

2500  RETURN  2  ' 

2510  FNENO 

2520  !  -  -•  - 

2530  DEF  FNP(Xin): 

2540  COM  /Consts/  Pl,Bl*) 

2550  IF  Xin>0  THEN 
2560  X-Xin 

2570  ELSE 
2580  X— Xin 

2590  END  IF 
2600  T-1/(1*P1*X) 

2610  Temp-0 

2620  FOR  1-5  TO  0  STEP  -1 
2630  Temp»T*Temp^B( I ) 

2640  NEXT  I 
2650  P-l-FN2CX)*Temp 
2660  IF  Xin<0  THEN  P-l-P 

2670  RETURN  P 

2680  FNENO 


